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^ ' Abstract 

^ ■ 

r^ ' In ultra-relativistic heavy ion collisions, the matter formed shortly after the collision is 

^v^ , a dense, out of equilibrium, system of gluons characterized by a semi-hard momentum scale 

1/^ ' Qs- Simple power counting arguments indicate that this system is over-occupied: the gluon 

occupation number is parametrically large when compared to a system in thermal equilibrium 

with the same energy density. On short time scales, soft elastic scatterings tend to drive the 

system towards the formation of a Bose-Einstein condensate that contains a large fraction of 

the gluons while contributing little to the energy density. The lifetime and existence of this 

condensate depends on whether inelastic processes, that occur on the same time scale as the 

elastic ones, preferably increase or decrease the number of gluons. During this overpopulated 

K^ , stage, and all the way to thermalization, the system behaves as a strongly interacting fluid, even 

» I ' though the elementary coupling constant is small. We argue that while complete isotropization 

C^ ' rnay never be reached, the system may yet evolve for a long time with a fixed anisotropy 

between average longitudinal and transverse momenta. 

1 Introduction 

One of the central theoretical issues in the description of heavy ion collisions is to understand 
how the partons that are freed by the collisions evolve into a thermalized system amenable to an 
hydrodynamical description. Let us recall that most of the produced partons originate from the 
small X components of the wavefunctions, that are dominated by gluon saturation and occupation 
numbers of order l/ag [1~3]. Such wavefunctions are well described by the Color Glass Condensate 
(CGC) effective field theory [4]. This effective theory allows in particular for the calculation of 
the energy-momentum tensor immediately after the collision. Because the chromo-electric and 



chromo-magnetic fields immediately after the collision are collinear to the collision axis, a config- 
uration called "Glasma" [5], the energy-momentum tensor is of the form T'^'^ = diag (e, e, e, — e), 
and therefore has a negative longitudinal pressure at very early times. Such an anisotropy be- 
tween the transverse and longitudinal pressures precludes a direct use of this energy-momentum 
tensor as initial condition for hydrodynamics, because the matching between the Glasma and the 
hydrodynamical evolution would require viscous corrections that are as large as the ideal terms. 

However this particular form of the energy momentum tensor, and the underlying structure 
of the fields, is not expected to last for a period of time much longer than 1/Qs, where Qs is 
the saturation scale. In fact, instabilities of various kinds [6] may lead rapidly (over the same 
time scale l/Qs) towards an isotropic energy momentum tensor. But even if isotropization of 
the energy-momentum tensor indeed occurs one outstanding issue remains, namely whether the 
phase space distribution functions relaxes towards the equilibrium Bose-Einstein distributions or 
not. The "bottom-up" scenario [7] provides a systematic way by which this relaxation occurs as 
a result of hard elastic and inelastic collisions. In ref. [8], it was demonstrated that anisotropy 
driven instabilities can significantly alter this picture. In this paper, we would like to address this 
question from a different point of view. 

Our paper is motivated by the basic observation that initially the gluon density in the Glasma is 
parametrically large compared to the value it should have in a system in thermal equilibrium with 
the same energy density. In systems where collisions conserve the particle number, any over (or 
under) -population in the initial condition can be accommodated by the appearance of a chemical 
potential in the equilibrium distribution. However, as we shall show, the initial over-population 
in the Glasma is so large, when the coupling constant is small, that the maximum allowed value 
for the chemical potential, that is fi — m (where ?Ti is a medium generated mass for gluons) 
is insufficient to account for the excess of gluons. This tension may be resolved in part by the 
dynamical generation of a Bosc-Einstcin condensate, corresponding to a large occupation of the 
zero momentum mode, and in part by inelastic processes that in the long run tend to tame the 
particle excess. 

After a short time scale of order 1/Qs, at which instabilities in the Glasma are expected to 
isotropize the system, it should be possible to describe the Glasma with color singlet distribution 
functions for both the particle content and the condensate. We will consider a generic form for 
such a distribution that has the following features: (i) it is dilute above some hard scale A, (ii) 
the gluon occupation number is saturated at a value l/ag below some "coherence scale" Ag, and 
(iii) between these two scales the gluon distribution is inversely proportional to the energy. At 
early times. As ~ A ^ Qg. In the expanding system, both A and Ag decrease in time; Ag decreases 
faster than A, so that thcrmalization is achieved by depletion of high energy modes. We will 
estimate the time it takes for Ag to become of order agA, and will argue further that this coincides 
with the thcrmalization time: in this picture, thcrmalization corresponds to the situation where the 
distribution function has acquired the equilibrium shape for most of the relevant momenta, thereby 
maximizing the entropy. We will argue that during the evolution to thermalization, because of the 
high occupancy of the low momentum modes, the system remains strongly interacting, although 
the coupling constant is small. It may also develop a fixed anisotropy ratio of (pi) relative to /p^ V 



2 The overpopulated quark-gluon plasma 

In the CGC description of heavy ion colhsions, the gluons that contribute dominantly to the energy 
density are freed over a time scale of order tq ^ Q^^ , with Q^ the saturation scale [9]. These gluons 
have typical transverse momenta of order Qs , and an energy density 

eo-6(r^07i)^^. (1) 

The number of gluons produced per unit volume is given by 

no = n(T - Q;') ^ ^ , (2) 

so that the average energy per gluon is indeed sq/hq ^ Qs- In fact Qs is initially the only scale 
characterizing the system, and the initial phase distribution function, a dimensionless object, is of 
the form fo{p/Qs,xQs,tQs ~ 1). 

3/4 

One may characterize the initial distribution of gluons by the dimensionless combination uq Cq 
l/as ■ In comparison, in an equilibrated system of gluons at temperature T, 



T^ , ncq-T'. (3) 



and 



ncqCcq 



f /4 ^ 1 . (4) 



There is therefore a mismatch, by a large factor as ^ I (in weak coupling asymptotics, ag <C 1), 
between the value of ne~^' ^ in the initial condition and that in an equilibrated system of gluons. 
We interpret this mismatch as an "overpopulation" of the initial distribution. Since the gluons are 
bosons, it is natural to explore the possibility that the system copes with this overpopulation with 
a new equilibrium state involving a Bose condensate. 

At this time, we assume that the equilibrium state is reached by elastic processes that con- 
serve the number of gluons, and consequently introduce a chemical potential in the equilibrium 
distribution function that one is looking for. The energy density and number density read then 

/cq(p) , "cq = / /cq(p) • (5) 



cq I LOp iQqyfij , ;tcq — / J cq 

' p J p 

where 

The temperature T — 1//3 and the chemical potential are adjusted so as to reproduce the initial 
values of eo and no- Here Wp is the energy of a gluon with momentum p. An important feature 
of any dense system of gluons is that, as a result of their many-body interactions, gluons develop 
effective medium dependent masses, that is, ujp^Q = m, 7^ . For instance, in a weakly interacting 

system of gluons in thermal equilibrium, this mass can be obtained in the well known Hard Thermal 

1/2 
Loop (HTL) approximation [10], where m ^ as T. For the initial distribution, we estimate the 

gluon mass as 



whereas, when the plasma has reached equiUbriuni, m ^ as T ^ a J Qg, with T ^ eoq ~ 
Qs/da ■ Note that in this discussion we shah make no distinction between the mass defined from 
the spectrum, m = Wp^oj and the screening mass. Both are parametrically comparable, and we 
shall often refer to m as the "Debye mass" . 

Since /oq(fc) is a growing function of the chemical potential, one way to cope with an excess of 
particles is to have a positive chemical potential. Note however that the chemical potential cannot 
grow larger than to, or else /oq would become negative. There is therefore a maximum number 
density that can be accommodated by the introduction of a chemical potential at a given T, 

d^k 1 ^3. (3) 



(27r)3 e/^("fc-'") - 1 

In this last estimate, we have used that to ^ T, so that the integrals are dominated by the region 
^^ ^ P ^ T where T/{ujp — to) ~ T/p. This calculation reveals that rimax ^ Qt/cts remains 
parametrically smaller than the initial gluon density no ^ Q^/ag. 

We are then led to the conclusion that, when the gluons undergo only elastic collisions, Bose 
condensation occurs, with the equilibrium distribution function taking the form 

/eq(fc)=^c^(fc)+ ^,(,J„„)_^ , (9) 

with ric the density of particles in the condensate, defined as the difference between the total 
density and the density of "thermal particles" (i.e., rimax). In fact, most of the particles are to be 
found in the condensate. We have indeed 

ne^^(l-ay^). (10) 

Note however that the condensate particles carry only a small fraction of the total energy density, 
since the energy of the condensed particles is 

n.m^^al^^Qs^al/^T^^eo. (U) 

The formation of a condensate is intimately associated with particle number conservation. 
When inelastic processes occur sufficiently rapidly, the number of gluons is not conserved anymore. 
Under such circumstance, no chemical potential can appear in the distribution function, and neither 
can a singular component representing the Bose condensate develop. The only equilibrium function 
is of the form of Eq. (6) with ^ = 0. In order to reach this equilibrium distribution, the system 
has to decrease its number of particles, via inelastic processes. 

It is also instructive to look at the entropy of the system, s ^ J In /p. This is dominated by 
the hard momenta, so that initially, s ^ Q^ (to within logarithmic corrections). In the case where 
there is no condensate, the equilibrium entropy is simply s ^ T^ ^ Q^/as . As expected, there is 
entropy increase as the shape of the momentum distribution evolves from the initial distribution 
towards the thermal distribution. Note that this increase is accompanied by the decrease of the 
particle number, so that in equilibrium s ^ T^ ^ n. When elastic collisions dominate, a condensate 
forms. The condensate carries no entropy, and absorbs the excess particles. The equilibrium state 
is characterized by the same equilibrium entropy s ^ T^, but now this entropy is carried only by the 



thermal particles, that is, s ^ nmax- In both cases, when equilibrium is reached the overpopulation 
disappears, as it should, namely, rige"^'* ^ 1, with rig ^ n when there is no condensate and 
Tig ^ rimax in the presence of a condensate. 

The thermodynamical considerations of the present section lead us to expect two possible 
equilibrium states, given the initial condition. Either a system with a Bose condensate, if the 
approach to equilibrium is driven by elastic collisions, or a system with fewer number of particles 
if inelastic processes are important. Note however that the presence of inelastic, particle number 
changing, processes does not preclude the possibility that a transient condensate develops during 
the evolution of the system. This is a dynamical issue that depends on the respective rates of 
particle production versus particle annihilation processes. We explore this question in the next 
sections. 

3 Kinetic evolution dominated by elastic collisions 

In order to address the question of how the system evolves towards its equilibrium state, we shall 
rely on a simple kinetic description based on the following transport equation 

dtf{k,X)^Ck[f], (12) 

where Ck[f] is the usual collision integral. We ignore at this point all drift terms in the left 
hand side. We assume that initially the system is isotropic, a property that is preserved by the 
evolution. In this section we focus on a gluon system in a non-expanding box. The effect of 
longitudinal expansion will be discussed later. We also assume that gluons undergo only elastic 
collisions, deferring the discussion of the effect of inelastic number changing processes to the next 
section. 

Our main goal is to understand how collisions drive the initial distribution towards local equilib- 
rium, and get a measure of the basic time scales involved. A detailed answer can only be obtained 
through explicit numerical solutions of the Boltzmann equation. Results of such calculations will 
indeed be presented elsewhere [11]. Here, we shall argue that we can capture the dominant qualita- 
tive features of the solution by assuming that the evolution is dominated by only two scales, Ag and 
A. We shall assume that the 1/uJp thermal distribution gradually builds up from energy Ag where 
the distribution function is ~ l/ttg down to A. To be concrete, although such an explicit form is 
not really needed in our arguments, we may assume that at all times t > 1/Qs, the distribution 
function takes the form 

f{p) for p < As, fip) ^ for As < p < A, f{p) - for A < p. (13) 

as as Wp 

At i ~ l/Qs, both scales Ag and A coincide with Qs- As time progresses, the two scale separates, 
with As decreasing quickly, and A evolving much more slowly. Thermalization is reached when 
Ag/A ^ ttg, at which point, /(A) becomes of order unity. 

A more precise definition of these two scales can be obtained by looking more closely at the 
collision integral. In the small angle approximation, assuming 2 — )■ 2 elastic scattering and isotropy, 
standard manipulations lead to [11] 



dl 
dt 



coll P 



l + l^,/(p)(l+,/(rt) 



(14) 



The fixed point solution of this equation is a Bose-Einstein distribution with temperature T = 
As/as (and indeed at thermahzation, T ^ A ~ Ag/ag). The two scales Ag and A may be obtained 
from the integrals 

as Jo dp ai Jq 

Remarkably, in the regime where f ^ ^ {f ^ l/o^s), all dependence on as drops from the collision 
integral. 

By taking moments of the collision integral above with arbitrary powers of p, it is not difficult 
to show that the typical collision time is given by 

A 

tscat — -T^, (16) 

which is independent of ag. This collision time should not be confused with the thermahzation 
time that we shall define later, and which depends on as- The collision time tgcat is also a function 
of time that we shall determine shortly. 

In all parametric estimates to be done below, we shall exploit the fact that the integrals are 
dominated by the largest momenta, of order A, large compared to the Debye mass, so that the 
distribution in the interesting region can just be taken to be f{p) ^ As/{asp) up to a cut-off of 
order A. Thus the number of gluons associated with this distribution is simply 

ng - — A^A,. (17) 

as 

In addition, one has the contribution of the condensate, ric, 

n = nc + ng, (18) 

with initially n = hq and ric — 0. Similarly, the energy density in gluon modes is 

, ^ IAsA^, (19) 

as 

to which should be added the energy of the condensate 

ec ^ m Tic (20) 

in order to get the full energy density. To determine Ec we need to estimate m, 

m'^as I dpp'^ ^AAs, (21) 

so that 

cc ~ nc^/AJi. (22) 

Initially, when Uc is small, ec represents a small correction to the energy density carried by the 
gluons. As times goes on, ric increases but Cc remains small, as we shall verify. 



and 



Our determination of the time dependenee of the scales Ag and A wiU rely on energy conserva- 
tion, 

AsA^ -- constant . (23) 

as well as the simple estimate for the scattering time given above, Eq. (16). The scattering time is 
itself a function of time. It is natural in the present context to look for a power law dependence, 
Qstscat ^ (Qsi)"- Then a simple analysis of the moments of the kinetic equation (12) reveals that 
the only sensible choice is a = 1 (provided one is not too close to equilibrium). We therefore set 

tscat ~ t. (24) 

With this assumption about the time dependence of the collision time, and imposing energy 
conservation, one easily determines the evolution of the two scales. We get 

3 

As - Qs (j) ' (25) 

A ^ Q. (^) ' (^^^ 

The number density of gluons rig ^ A^A,, decreases as ^ (to/i)^^^, while the energy carried by 
the gluons, ^ AgA^ remains approximately constant. The Debye mass decreases slowly in time, 
m ^ Qsito/ty^''' , so that indeed the energy carried by the condensate particles, with density 
"n-c ~ no[l — (to/i)^^'^], remains negligible. 

The thermalization time, determined from Ag ^ agA, is 

ith-^f-)'. (28) 

Qs \asj 

We notice that since the scale A is increasing in time, and therefore so is the entropy density s ~ A'^, 
there is entropy generation during the evolution, s ^ A'^ ^ Qs (V^o)^^^- When t = ith, s ~ Q^/as , 
which is the equilibrium entropy ^ T'^. Note also that our expression (24) for the scattering time 
interpolates between 1/Qs, the scattering time in the initial plasma, and Os Q^^ ^ ^/{'^s'^) in 
the equilibrated plasma. In fact, near thermalization, the scattering time may be given a familiar 
interpretation in terms of kinetic theory: 1/iscat ^ o'n, where a ^ a^/A^, and n ^ A'^, which yields 
indeed 1/tscat ~ ci^s^- Such an interpretation does not hold in the initial state where 1/tscat ^ Qs, 
with no reference to the coupling constant. 

Finally, the effect of quarks can be estimated. Quarks should have a phase space density of 
order 1 up to the scale A. Therefore the number density of quarks is 

^quarks ^^ ^^ \^^ ) 

Note that initially the quark number density is one order of as smaller as compared with that of 
the gluons rig ~ AgA^/ag. But by the thermalization time when Ag ~ agA, they are of the same 
order and become equally important only at this time. Quarks may then, to first approximation, 
be ignored until the time of thermalization. 



4 The effects of inelastic processes 

Inelastic particle production or annihilation processes will modify the collision integral on the right 
hand side of the transport equation (12). Consider for example, the contribution of an n — > m 
process to the collision term. The vertices contribute a factor a"+'"~^. There is a factor of 
(As/q!s)"^™^^ arising from the distribution functions (one factor for each distribution, except the 
one whose momentum one is following; besides, the products containing n + m factors / cancel 
between the gain and loss terms). It follows that the coupling constant disappears explicitly. 
Furthermore, as shown for instance by Mueller et al. [12], there is an overall infrared singularity in 
the multiparticle production diagrams, which is cutoff by the Debye scale, and this yields a factor 
of (1/to^)""'"'""'*. Using m^ ~ AgA, we obtain overall a factor of Ag/A"+™~'*. This is balanced by 
a factor emerging from the remaining phase space integral. Since this integral is infrared finite, it is 
dominated by momenta of order A and is therefore proportional to a power of A. For dimensional 
reason, this is A"''"™"^, leaving an expression for the inelastic scattering which is parametrically 
identical to that for elastic scattering^, namely tgcat ~ ^/^s- 

Our procedure for determining Ag and A above used only igcat ^ t and energy conservation. 
Therefore, quite remarkably, including the effects of inelastic scattering does not change the scaling 
behaviour for Ag and A. 

There are modifications of the treatment of the condensate however. As we have already 
mentioned, inelastic processes will inevitably lead to an equilibrium state without a condensate. 
The question then arises of whether such a condensate can exist as a transient state for a sufficient 
amount of time to infiuence the dynamics of the system. Of course, the answer to this question can 
only be obtained after a detailed numerical analysis of the solution of the transport equation. We 
can however offer the following lines of reasoning. In the small angle approximation to the transport 
equation for elastic processes, one finds, in case of overpopulation, that the gluon distribution 
develops very rapidly a 1/p behavior near p = 0, which eventually generates the delta-function 
singularity characteristic of the condensate. Unless the inelastic processes change this singular 
behavior of the distribution at small p, the elastic contribution to the collision integral provides 
therefore a source term for the condensate. There may be other source terms associated with 
higher order, multi-particle, processes. Inelastic scattering terms will also contribute a sink term. 
Since elastic and inelastic processes evolve with the same time scale, it is conceivable that a balance 
between the source and sink term can be achieved, and that a condensate is created and survives 
for most of the evolution till thcrmalization. 

5 Effect of the longitudinal expansion 

An important feature of the matter produced in ultra-rclativistic heavy ion collisions is its strong 
longitudinal expansion. Assuming longitudinal boost invariance, and focussing on the central slice 
z — 0, one may capture the main effect of this expansion by adding to the left hand side of the 
kinetic equation a drift term of the from 



t ^' dt 



= C[f] , (30) 



^We thank Guy Moore for pointing out a mistake in this analysis in an earlier version of this manuscript. 



where the notation -^ stands for a time derivative at constant p^t. In the absence of the colhsion 

term, this equation admits free streaming solutions of the form f{p±,Pzit) — f{pj_,pzt/to)- Thus, 
as an immediate effect of the longitudinal expansion, an isotropic initial distribution will flatten 
in the z-direction, on a time scale of order tg ~ l/Qs- This is an effect that will potentially delay 
complete three dimensional isotropization of the particles. We shall come back to this issue in the 
next section. Here, we shall make a simplifying assumption that the system may evolve under the 
combined effect of longitudinal expansion and collisions with a fixed anisotropy. We shall quantify 
shortly the degree of anisotropy. 

By integrating over momentum the kinetic equation multiplied by the energy one obtains 

dte + '-^ = , (31) 

where e is the energy density, £ — J (2-n{i ^pfpi ^^'^ Pl the longitudinal pressure, Pj^ — J i^tAa zf'fp' 
One may analyze the effects of longitudinal expansion by parameterizing the longitudinal pressure 
in terms of energy density, namely, by assuming P^ — 5 e where the multiplicative factor 5 can be 
in the range [0, 1/3] with (5 = corresponding to the completely free streaming case and 5=1/3 
corresponding to ideal hydrodynamic expansion after isotropization. Of course, the assumption 
that 5 is independent of time is a strong assumption. We make it here in order to focus on the issue 
of how collisions redistribute momenta and thereby generate the shape of a thermal distribution 
for an expanding system. Note also that the nature of the local equilibrium changes somewhat as 
a function of 5. For instance, when (5 = the local equilibrium is essentially two-dimensional, a 
situation that precludes the formation of a Bose condensate. 

Within the present assumption the equation (31) becomes an equation for the evolution of the 
energy density 

^.(i) - <to) (y) . (32) 

This equation, together with our previous estimate of the collision time (we can verify that the 
linear relation tscat ~ t remains unaffected by the expansion), yields the following estimates 

. X (4+5)/7 / X (l+25)/7 



From these, one easily obtain the estimates of the gluon density, and of the Debye mass 

The thermalization time, obtained as before from the condition Ag = agA, is given by 

'tth\ /I 



to J V tts / 

By comparing to the static case studied previously, we see that the expansion has the effect of 
delaying thermalization slightly. (Formally, one may recover the static case by setting 5 = —1, 
which corresponds to constant energy density.) 



At this point, we may consider two scenarios. First, we assume particle number conservation. 
Then by integrating the transport equation over momenta, and noticing that the colhsion term 
does not contribute, we get 

atn + ^=0, n = nof|V^i. (36) 

t \ t / as t 

Note that the effect of the expansion is to decrease the parameter ne^^^^ that characterizes the over- 
population (ignoring thermalization processes beyond those responsible for maintaining isotropiza- 
tion): 

-3/4 / '^O \ I to 



Clearly, for S = 1/3, which corresponds to isotropic expansion, the decrease of the density and 
that of the energy density combine so as to leave the overpopulation parameter unchanged, while 
a fast decrease is achieved for the free streaming case 5 = 0. For moderate values of the anisotropy 
(more precisely for 6 > 1/5), a condensate can form, with density 



^ (-l+5<5)/7- 



(38) 



As before, we can verify that the energy carried by the particles in the condensate is subleading: 

(39) 



., ^(5-115)/14 



Eg V t 

which decreases with increasing t. 

Alternatively, we may assume that particle number is not conserved. The analysis of this case 
follows that of the non expanding box, with the same uncertainty concerning the final role of 
inelastic scattering. 

6 The asymmetry 

At early times, a mechanism that attempts to restore isotropy between p^ and p^ are the Weibel 
instabilities [15-18,8,6] .^ These arise from purely imaginary screening masses in the presence of an 
asymmetric momentum distribution of gluons [19,20]. Generically, the time scale for the restoration 
of such symmetry is the inverse of the Debye mass, which we have argued to be m ~ y/AA^. Note 
that 1/m ^ l/VAAg ^ A/Ag = iscat for times when A ^ Ag. This of course confirms that the 
time to restore isotropy is much less than the scattering time. If isotropy were maintained, one 
would require that 6 = ^; the system would then evolve in time according to ideal hydrodynamics. 
However, as the system evolves in time, the Weibel instability, that operates on soft modes with 
momentum p < to is likely not efficient enough at isotropizing hard modes with p ^ A ^ to,. It is 



^We will distinguish these instabilities from the initial "leading" instabilities that were responsible for isotropizing 
the system in the first place on a time scale of order 1/Qs [14]. We however leave open the possibility that the 
subsequent dynamical evolution of the quantum system in ref. [14] may overlap at later times as well with the 
physics under discussion here. 
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possible that there may be some other mechanism similar to the Weibel instability that operates 
on higher momentum scale and generates isotropy there. 

We may argue that scattering, while probably not sufficient to fully restore isotropy, may 
nevertheless maintain the system in a state of fixed anisotropy for a long time. For the sake of 
illustration, we can write down a simple kinetic equation that achieves this goal. To this aim, let 
us integrate the left hand side of the kinetic equation after multiplying it by either pi or p\. One 
gets 



2 df 



.2 ^ 

dt 



= dM 



1 



t 



(pi), 



(40) 



= dt{pi) + -{pi), I VI 

with (jij,) = J pI f, and similarly for {p]_)- Defining the asymmetry p^ in the momentum distri- 
bution by 

pI^pI-IpI, (41) 

with pI — {p^ + 2p^)/3, and using the equations above for {p^) and (p^), we get 



2 df 



Pzt 



-0,{pl) + l-^{pl) + l-^{p^). 



(42) 



It is plausible to assume that the collision term will contribute a "relaxation" force for the asym- 
metry. We therefore complete the equation as follows 



dt{pl) + ^yj + l^{p') = -^{pl), 

with K a constant characterizing the strength of the collisions. 
We will look for scaling solutions such that 



{?') - {p\ 



toY 



(43) 



(44) 



and further introduce a time-dependent dimcnsionless parametrization of the asymmetry /p^ ^ , 

{pl) = m{p'). (45) 

We find 

tdtC+(^ + K-7i\c+^^0. (46) 

As anticipated, this equation has a solution which relaxes towards a fixed value of C, namely 

2/3 



'' 7/3 + K-^' 

One may relate (approximately) C to the parameter S introduced earlier: 

S^{pl)/{p') = il + 20/S. 
One can then eliminate 77 = (8 + 9S)/7 and get C as a function of k: 

16 -f 21k 



c 



36 



1008 
(16 + 21k) 



(47) 



(48) 



(49) 
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Depending on the strength of the colhsions, represented here by the parameter k, various (negative, 
i.e. (pI) < (p5_)) values of the anisotropy can be reached, from C sa —1/2 for smaU k (in fact we 
must keep k > 1/7 for S to stay positive) to C = when k — > cx). 

7 Suminary 

This paper argues that the Glasma formed in the early stages of heavy ion collisions is strongly 
interacting with itself up to parametrically late times when the system thermalizes. In particular, 
we show there are scaling solutions to the transport equations from which the coupling constant 
has disappeared. In addition, there may exist a transient component of the system, which is a 
Bose-Einstein condensate. If this scenario is realized, it may have a profound impact on the way 
in which we describe the properties of the Quark-Gluon Plasma in heavy ion collisions. 
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